Astronomy & Astrophysics manuscript no. 0328 


February 2, 2008 


(DOI: will be inserted by hand later) 





Line formation in solar granulation 

IV. [Oi], Oi and OH lines and the photospheric O abundance 

Martin Asplund^, Nicolas Grevesse^'^, A. Jacques Sauval^, Carlos Allende Prieto^ and Dan Kiselman^ 

, ^ Research School of Astronomy and Astrophysics, Mt. Stromlo Observatory, Cotter Rd., Weston, ACT 2611, Australia 

' ^ Centre Spatial de Liege, Universite de Liege, avenue Pre Aily, B-4031 Angleur-Liege, Belgium 
^SJ , Institut d'Astrophysique et de Geophysique, Universite de Liege, Allee du 6 aout, 17, BSC, B-4000 Liege, Belgium 

^ . * Observatoire Royal de Belgique, avenue circulaire, 3, B-1 180 Bruxelles, Belgium 

^ ' ^ McDonald Observatory and Department of Astronomy, University of Texas, Austin, TX 78712-1083, USA 

^\ , ^ The Institute for Solar Physics of the Royal Swedish Academy of Sciences, AlbaNova University Centre, SE-106 91 

■ Stockholm, Sweden 

Accepted for Astronomy & Astrophysics 

^ , Abstract. The solar photospheric oxygen abundance has been determined from [O I], O I, OH vibration-rotation and OH pure 

• rotation lines by means of a realistic time-dependent, 3D, hydrodynamical model of the solar atmosphere. In the case of the 

I O I lines, 3D non-LTE calculations have been performed, revealing significant departures from LTE as a result of photon losses 

■ in the lines. We derive a solar oxygen abundance of logeo = 8.66 ± 0.05. All oxygen diagnostics yield highly consistent 
CNj ' abundances, in sharp contrast with the results of classical ID model atmospheres. This low value is in good agreement with 

measurements of the local interstellar medium and nearby B stars. This low abundance is also supported by the excellent 
correspondence between lines of very different line formation sensitivities, and between the observed and predicted line shapes 
and center-to-limb variations. Together with the corresponding down-ward revisions of the solar carbon, nitrogen and neon 
, abundances, the resulting significant decrease in solar metal mass fraction to Z = 0.0126 can, however, potentially spoil the 

^ l' impressive agreement between predicted and observed sound speed in the solar interior determined from helioseismology. 
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1 . Introduction tain consistency is very restricted yet the discrepancy is present 

in full. 

Oxygen is the most abundant element in the Universe with a ^^^^ ^^^^^^^^ commonly adopted solar oxygen abun- 

H non-Big Bang nucleosynthesis origin. As a consequence, oxy- ^^^^^ ^ g gg ^ q ^^^^^^^^ ^ q^^^^^^^ ^^g^^ 

. - gen plays a central role in many different fields of astrophysics -ru- u- ^ ■ ^^ u- u u /i .ju i t 

° ^ ^ ^ ^ This historically high abundance was suggested by analyses of 

ranging from supernova physics and galaxy evolution to dat- ^j^^ forbidden [Ol] 630.0 nm line (Lambert 1978) as well as 

ing stars and production of the light elements through cosmic r\u u » a » i- ,.u ■ t a 

^ f fe & OH vibration-rotation and pure rotation lines in the infrared 

ray spallation. Yet it appears that in many crucial objects for (o^eyesse et al. 1984; Sauval et al. 1984) using the ID hy- 

which accurate knowledge of the oxygen abundances is neces- ^^^^^^^.^ Holweger-Miiller (1974) semi-empirical model of the 

sary the oxygen content is hotly debated. Recent disputes re- , . ■ j t t-c r f /~v »u »u u j 

^ ^ solar atmosphere and LTE line formation. On the other hand, 

volve around the overabundance of oxygen in metal-poor halo uj • •j ^juti, -^^juu 

^ a much lower abundance is indicated by the permitted high- 
stars (see Asplund & Garcia Perez 2001; Nissen et al. 2002 .... r »u <.u m <. • i <. <. ^^^ 

^ ^ excitation O I lines, most noteworthy the IR triplet at 777 nm, 

and references therein), the Galactic radial abundance gradient ■ t ■ <.u j i » u -.u n-tr 

' ° when employing the same model atmosphere with non-LTE 

(Rolleston et al. 2000; Cunha & Daflon 2003), and, astonish- ,. f .■ r^u- a- * u » m j u » A t 

^ /' ' jjjjg formation. This discrepancy of about 0.2 dex between dif- 

ingly, the solar oxygen abundance. Partly these disagreements f »uj ■ a- , u f»u ui a 

^ •'^ ^ ferent abundance indicators have often been blamed on over- 
stem from differences in the adopted input data (e.g. Q f -values, ^- ^ A A » t 1 1 »U A ■ -lU • 

^ ^ V fe i/j ' estimated departures from local thermodynamic equilibrium 

effective temperatures T.ff, surface gravities log 5) but more ^ltE) for the Ol Unes (Tomkin et al. 1992; Takeda 1994). 

importantly they reflect the choice of spectral lines to derive j ■ j -.u rT-c u a f»u » • 1 » • 1 * »u r/^xi a 

J J f Indeed, the LTE abundance of the triplet is close to the [O I] and 

the abundances using classical ID stellar model atmospheres, u j u j -^u »u o 1 N/t- n j 1 c u 

° ^ OH-based abundance with the Holweger-MuUer model. Such 

In particular in the solar case, the freedom of parameters to ob- 



On the customary abundance scale defined as e(X) — 10^ 
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argumentation, however, ignores the available observational 
evidence such as the center-to-limb variation that the triplet is 
definitely not formed in LTE (Altrock 1968; Sedlmayr 1974; 
Kiselman 1993; Kiselman & Nordlund 1995). Uncertainties 
surround not only the O I lines but also the other indicators. The 
[O l] line is very weak and blended with a Ni I line (Lambert 
1978; Reetz 1998; Allende Prieto et al. 2001). The OH lines 
are, like all molecules in the relevant temperature regime of the 
solar atmosphere, very temperature sensitive and therefore sus- 
ceptible to surface inhomogeneities like granulation (Kiselman 
& Nordlund 1995; Asplund & Garcia Perez 2001), which are 
ignored in classical ID model atmospheres. Unfortunately, no 
guidance regarding the solar oxygen abundance is available 
from meteorites, contrary to the situation for most other ele- 
ments (Grevesse & Sauval 1998). Like carbon and nitrogen, 
oxygen is a highly volatile element which has partly escaped 
from the most primitive meteorites, the Cl-chondrites, leaving 
only a fraction of the original amount (10%, 3% and 54% of C, 
N and O, respectively). 

Recently, Allende Prieto et al. (2001) suggested a signif- 
icant down-ward revision of the solar oxygen abundance to 
logeo = 8.69 ± 0.05 based on a study of the [Ol] 630.0 nm 
line using a time-dependent 3D hydrodynamical model of the 
solar atmosphere (Asplund et al. 2000b). Most of the 0.24 dex 
difference with the previous higher value is attributed to a pre- 
viously ignored Ni I blend (—0.13 dex), although the difference 
between 3D and ID models compounds the effect (—0.08 dex). 
In addition, a revision of the transition probability of the line 
causes a further 0.03 dex decrease. A similar low value was 
obtained by Holweger (2001) driven primarily by permitted 
Ol lines: logeo — 8.74 ± 0.08. This new low O abundance 
largely resolves another puzzling problem with the old high 
oxygen abundance, namely that the Sun is apparently oxygen- 
rich compared with the present-day local interstellar medium 
(e.g. Meyer et al. 1998; Andre et al. 2003) and nearby B stars 
(e.g. Cunha & Lambert 1994; Kilian et al. 1994; Sofia & Meyer 
2001). Since the birth of the Sun some 4.5 Gyr ago, the gas 
in the solar neighborhood from which these hot stars have re- 
cently formed should have been further chemically enriched by 
nuclear-processed ejecta from dying stars. Several explanations 
like migration of the solar Galactic orbit have been proposed to 
address this conundrum but they may all be superfluous if the 
new low solar oxygen abundance is confirmed. 

Relying solely on a single line for an as important quantity 
as the solar O abundance is no doubt highly unsatisfactory, in 
particular given the significant Ni blend of the [O l] 630.0 nm 
line. It is therefore clearly of importance to extend the work 
by Allende Prieto et al. (2001) to include also the O I and OH 
lines analysed with the same 3D hydrodynamical model atmo- 
sphere. Such a study is presented here, which firmly establishes 
the reality of the low solar photospheric oxygen abundance. 
It is noteworthy that for the first time ever the different line 
indicators give consistent results, which strongly supports the 
high degree of realism obtained with the new generation of 3D 
model atmospheres suggested by previous studies. 



2. Analysis 

2. 1. Atomic and molecular data 

[O I] lines: The g/-value of the [O l] 630.03 nm line has re- 
cently been revised to log 5/ = —9.717 (Storey & Zeippen 
2000), which we adopt here. This is a slight revision from 
the previously commonly adopted value of log gf = —9.75 
(Lambert 1978). The 630.03 nm feature has a non-negligible 
contribution from a Nil line (Lambert 1968; Reetz 1998; 
Allende Prieto et al. 2001) which must be taken into account. 
Johansson et al. (2003) has very recently measured the tran- 
sition probability of the Ni line and found loggf = —2.11. 
However, since only the product eNi • gf enters the x^-analysis 
of the [Ol] 630.0 nm line profile fitting as a free parameter 
(Allende Prieto et al. 2001), this new value will not modify 
the derived O abundance but only the solar Ni abundance sug- 
gested by this one line. The isotopic splitting for the Ni line also 
determined by Johansson et al. (2003) is taken into account for 
the synthesis of the 630.03 nm line. The corresponding data for 
the weaker [Ol] 636.37 nm line is log 5/ = —10.185 (Storey 
& Zeippen 2000). In the analysis of this transition care must be 
exercised due to problematic blending of CN lines and a Ca I 
auto-ionization line at the wavelength of the [Ol] 636.37 nm 
line. 

O I lines: As described below, we chose to retain only six 
permitted Ol lines (615.81, 777.19, 777.41, 777.53, 844.67 
and 926.6 nm) for the abundance analysis to minimise the er- 
rors stemming from blends, uncertain g/-values, continuum 
placement and too weak spectral lines. The transition proba- 
bilities and lower level excitation potentials for the five lines 
were taken from the NIST database^. The VALD database^ 
(Piskunov et al. 1995; Kupka et al. 1999) was consulted for ra- 
diative broadening and central wavelengths. Collisional broad- 
ening by H (commonly referred to as van der Waals broad- 
ening) was computed from the quantum mechanical theory 
of Anstee & O'Mara (1995), Barklem & O'Mara (1997) and 
Barklem et al. (1998). This removes the need for the con- 
ventional enhancement factors to the Unsold (1955) classical 
recipe for collisional broadening. Some of the relevant line 
properties are listed in Tabled 

OH vibration-rotation and pure rotation lines: We selected 
the best 70 vibration-rotation lines of the (1,0), (2,1) and (3,2) 
bands as well as the best 127 pure rotation lines belonging to 
the 0-0, 1-1, 2-2 and 3-3 bands in the infrared solar spectrum. 
The dissociation energy of OH is very accurately known since 
quite a long time: Dq (0H)=4.392 ± 0.005 eV (Carlone & 
Dalby 1969). The molecular partition functions and equilib- 
rium constants have been taken from Sauval & Tatum (1984). 
Lower level excitation potentials were adopted from Melen et 
al. (1995). The ^/-values were calculated by E. van Dishoeck 
(private communication) adopting the Electric Dipole Moment 
Function (EDMF) of Nelson et al. (1990). We checked that this 
EDMF leads to more consistent results than other calculated 
EDMFs. We note that this EDMF has also been selected by 
Goldman et al. (1998) in their recent database for the OH X 

^ http://physics.nist.gov/cgi-bin/AtData/ main_as5l 
^ _http://www.as tro.uu.se/~vald. 
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ground state. A direct comparison of the ^/-values with those 
of Goldman et al. (1998) reveal essentially identical results. 

2.2. Observational data 

For the analysis of the forbidden [O l] and permitted O I lines 
which are all located in the optical region (600-850 nm) the so- 
lar flux atlas of Kurucz et al. (1984) has been employed. This 
atlas has a signal-to-noise ratio in the optical of 200-3000, and 
a resolving power of about 500,000. Minor adjustments to the 
tabulated continuum level were made for the small wavelength 
regions directly surrounding the lines. For the 3D analysis, the 
oxygen abundances suggested by the individual lines were de- 
termined by profile fitting which is possible in view of the in 
general excellent agreement between predicted and observed 
line profiles in 3D (e.g. Asplund et al. 2000a,b,c). This is also 
the main reason why we employ flux profiles rather than the 
computationally less demanding disk-center intensity profiles 
for the oxygen lines: the agreement between predicted and ob- 
served profiles is slightly worse in intensity for the O I lines. 
Their very significant departures from LTE modify the inten- 
sity profiles more than the corresponding flux profiles. We em- 
phasize, however, that very similar oxygen abundances are ob- 
tained for the O I and [O l] lines with flux and disk-center in- 
tensity profiles. 

The OH vibration-rotation and pure rotation lines are lo- 
cated in the infrared. Our analysis has made use of the 
Spacelab-3 Atlas-3 ATMOs'* solar disk-center intensity IR ob- 
servations recorded in November 1994 which are of better qual- 
ity (i.e. less telluric absorption) than in the IR atlas (Farmer & 
Norton 1989; Farmer 1994). The relevant OH lines were iden- 
tified and their equivalent widths measured (in milli-Kayser = 
milli cm"^, subsequently converted to the corresponding val- 
ues in wavelength units [pm]). 

Solar observations of the center-to-limb variation of several 
spectral lines, were carried out in October 22, 1997, with the 
Gregory Coude Telescope (GCT) and its Czerny-Turner echelle 
spectrograph (Kneer et al. 1987) at the Observatorio del Teide 
(Tenerife, Spain). The observations of the Ol 777 nm triplet 
have an estimated resolving power i? = 86, 000. Scattered 
light was estimated to amount to about 6% and corrected 
by comparing the observations at the center of the disk with 
the Fourier Transform Spectrograph atlas by Brault & Neckel 
(1987). The observations were done for disk-positions /i = 
1.00,0.97,0.87,0.71,0.50 and 0.17. The observed center-to- 
limb variations of the equivalent widths for the O I 777 nm 
triplet are very similar to those measured by Miiller et al. 
(1968) and Altrock (1968) except at the lowest /^-point, where 
our data is less steep. More details of the observations and data 
reduction are provided in a separate paper (AUende Prieto et 
al., in preparation). 

2.3. 3D hydrodynamical solar model atmosphere 

We employ the same 3D, time-dependent, hydrodynamical 
model of the solar atmosphere, which has previously been ap- 

* |http://remus.jpl.nasa.gov/atmos/ftp.at3.sun.htmI| 
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Fig. 1. The temporally and spatially averaged temperature 
structure of the 3D solar surface convection simulation (solid 
line) used for the 3D spectral line formation. The spatial aver- 
aging has been performed over surfaces of the same continuum 
optical depths at 500 nm. Note that the actual 3D simulation 
extends to much greater optical depths than shown here. Also 
shown are the temperature structures for the ID Holweger- 
Miiller (1974) semi-empirical solar atmosphere (dashed line) 
and the ID MARCS (Asplund et al. 1997) theoretical solar at- 
mosphere (dashed-dotted line). 



plied to solar spectrum line formation for the present series 
of articles (Asplund 2000; Asplund et al. 2000b,c; Asplund et 
al. 2004b,c) and elsewhere (AUende Prieto et al. 2001, 2002a, 
b; Asplund 2004). This new generation of model atmospheres 
has proven to be highly realistic, successfully reproducing a 
wide range of observational constraints, such as spectral line 
shapes, shifts and asymmetries (Paper I; AUende Prieto et 
al. 2002a), helioseismology (Rosenthal et al. 1999; Stein & 
Nordlund 2001), flux distribution and limb-darkening (Asplund 
et al. 1999b) and granulation properties (Stein & Nordlund 
1998). It is clear that the new 3D solar model atmosphere rep- 
resents a significant improvement over existing ID hydrostatic 
LTE models. 

In short, the equations of mass, momentum and energy 
conservation together with the simultaneous solution of the 
3D radiative transfer equation along nine inclined directions 
have been solved on a Eulerian mesh with 200 x 200 x 82 grid- 
points: for the spectral line formation calculation this was sub- 
sequently interpolated to a 50x50x82 grid. The physical di- 
mension (6 Mm horizontally and 3.8 Mm vertically with about 
1 Mm located above T500 = 1) of the numerical grid were 
sufficiently large to cover > 10 granules simultaneously and 
extends to nearly adiabatic conditions in the bottom. Periodic 
horizontal and open vertical boundary conditions have been 
employed. The simulations make use of the state-of-the-art 
MHD equation-of-state (Mihalas et al. 1988) and comprehen- 
sive opacities, including line opacities (Gustafsson et al. 1975 
with subsequent updates; Kurucz 1993). The effects of line- 
blanketing is accounted for through opacity binning (Nordlund 
1982) which is a form of multi-group opacities. Some rele- 
vant properties of the temporally and spatially (over surfaces of 



4 



M. Asplund et al.: Solar line formation: IV. The photospheric O abundance 



Table 1. Temporally and spatially averaged atmospheric stratification of the 3D solar surface convection simulation used for the 
3D spectral line formation presented here: temperature, density, gas pressure, electron pressure and vertical velocity (positive 
velocities inward) as well as their rms-scatter The spatial averaging has been performed over surfaces of the same continuum 
optical depths at 500 nm. Note that the original 3D model extends significantly above and below the optical depths for which 
the mean structure has been interpolated to and tabulated. We emphasize that it is not possible to explain all 3D effects by a 
comparison of such an averaged 3D model with existing ID models due to the absence of atmospheric inhomogeneities. 
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same optical depths) averaged 3D solar simulation is listed in 
Table[2 limited to the optical depth range relevant for spectrum 
formation. The mean temperature structure < T(t5oo) >t is 
shown in Fig.^ We emphasize, however, that the atmospheric 
inhomogeneities are also important for the spectrum formation, 
which has been fully taken into account in the 3D line forma- 
tion. The reader is referred to Stein & Nordlund (1998) and 
Paper I for further details of the construction of the 3D solar 
model atmosphere. 

For comparison purposes, we have also performed identi- 
cal calculations employing two well-used ID hydrostatic mod- 
els of the solar atmosphere: the semi-empirical Holweger- 
IVIuller (1974) model and a theoretical, LTE, line-blanketed 
MARCS model (Asplund et al. 1997). Their temperature struc- 



tures T(t5oo) are compared with the mean 3D stratification in 

Fig.m 

2.4. 3D LTE spectral line formation 

Equipped with the 3D hydrodynamical solar model atmo- 
sphere, 3D spectral line formation calculations have been per- 
formed for the [Ol], Ol and OH lines. The radiative transfer 
equation has been solved under the assumption of LTE for the 
ionization and excitation with the source function being equal 
to the local Planck function {Si, — B,j). Comprehensive back- 
ground continuous opacities (Gustafsson et al. 1975; Asplund 
et al. 1997) and an equation-of-state which accounts for excita- 
tion, ionization and molecule formation of the most important 
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elements (Mihalas et al. 1988) have been used, with instanta- 
neous chemical equilibrium being assumed for the molecule 
formation. For the [O l] and O I lines flux profiles have been 
computed from a solution of the radiative transfer equation 
along 17 different inclined directions (A^^ — N^p = 4 plus 
the vertical) making use of the periodic horizontal boundaries. 
For the OH lines disk-center intensity profiles have been ob- 
tained. The flux profiles have been disk-integrated assuming a 
solar rotation velocity of 1.8 kms^^. All 3D LTE spectral line 
formation computations have been carried out for 100 snap- 
shots covering 50 min of solar time after which temporal and 
spatial averaging have been performed. For each line, a local 
continuum flux/intensity has been computed with the same pro- 
cedure as for the line profiles, enabling proper normalized pro- 
files. Each line has been calculated with three different abun- 
dances differing by 0.2 dex from which the profile with the cor- 
rect equivalent width (OH lines) or line shape ([Ol] and Ol 
lines) has been interpolated. Extensive testing has shown that 
the procedure yields correct abundances to within 0.01 dex. No 
micro- and macroturbulence enter the 3D calculations as the 
self-consistently computed convective velocities and temper- 
ature inhomogeneities induce Doppler shifts which give rise 
to excellent agreement between observed and predicted line 
shapes without such additional line broadening which is always 
needed in standard ID analyses (Asplund et al. 2000b; Asplund 
2000; Allende Prieto et al. 2002a). 

The corresponding ID calculations based on the MARCS 
and Holweger-Miiller model atmospheres have been performed 
using exactly the same codes and procedures as for the 3D 
case, with the exception that a microturbulence of ^tuib = 
l.Okms"^ has been adopted. To investigate the sensitivities 
of the lines to this parameter, all ID profiles have also been 
computed using fturb — l.Skms^^. 

2.5. 3D no n- LTE spectral line formation for Oi 

It is well-known that the high-excitation permitted Ol lines 
suffer from significant departures from LTE (Sedlmayr 1974; 
Kiselman 1991; Kiselman & Nordlund 1995; Kiselman 2001). 
It is therefore paramount to investigate these non-LTE ef- 
fects also in 3D and include non-LTE abundance corrections 
in the final analysis of the Ol lines. The 3D non-LTE line 
formation calculations have been performed with MULTI3D 
(Botnen 1997; Botnen & Cai'lsson 1999; Asplund et al. 2003a), 
which iteratively solves the rate equations simultaneously with 
the radiative transfer equation assuming statistical equilibrium 
(dNi{x , y , z) / dt = 0). Due to the very computationally de- 
manding nature of 3D non-LTE line formation calculations, 
compromises must be made in terms of the number of snap- 
shots for which the non-LTE solution is obtained. From the 
50 min time sequence for which the 3D LTE line formation cal- 
culations were performed, two typical snapshots were chosen 
which were sufficiently separated in time to be considered in- 
dependent. To improve the numerical accuracy in the non-LTE 
case, the 50x50x82 grid used for the 3D LTE spectrum syn- 
thesis were interpolated to a finer 50 x 50 x 100 grid extending 
down to at least logrsoo = 2 for all vertical columns. Test 
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Fig. 2. Grotrian term-diagram of the employed 23-level O 
model atom with 22 bound states of O I and the ground state 
of O II. The forbidden [O l] 630.0 and 636.3 nm lines originate 
from the 2p^ ground state, while the Ol 777 nm triplet is 
between the high-excitation 3s '^S" and 3p '^P levels. The 43 
bound-bound transitions are drawn. All O I levels are connected 
by bound-free transitions with the ground state of O II but are 
not shown for clarity 



calculations were also performed on a smaller 25x25x100 
grid, yielding very similar results. Tests with different vertical 
depth-scales yielded insignificant differences in the resulting 
non-LTE abundance corrections. 

The radiative transfer equation was solved for 24 outgoing 
inclined rays using a short characteristic technique and mak- 
ing use of the periodic horizontal boundary conditions. The 
line broadening caused by Doppler shifts introduced by the 
convective motions in the hydrodynamical model atmospheres 
was taken into account. No microturbulence or macroturbu- 
lence thus enter the 3D calculations. The background contin- 
uous opacities and equation-of-state were computed with the 
Uppsala opacity package (Gustafsson et al. 1975 and subse- 
quent updates). For completeness we note that this equation-of- 
state is not the same as used in the construction of the 3D model 
atmosphere but very similar in the relevant temperature regime. 
The effect of this difference will however be minute since both 
the LTE and non-LTE line profiles for the two selected snap- 
shots used for the derivation of the non-LTE abundance correc- 
tions are both based on the same input physics in all aspects, 
as described in detail below. This non-LTE abundance correc- 
tion is then applied to the LTE estimates based on the whole 
100 snapshot long temporal sequence, which indeed have the 
same equation-of-state as the 3D model atmosphere. Finally, 
we note that all four codes employed in this study (3D hydro- 
dynamics, 3D LTE line formation, MULTI for ID non-LTE line 
formation and MULTI3D for 3D non-LTE line formation) rely 
on the Uppsala opacity package for continuous opacities. 

Line-blanketing using data from Kurucz (1993) were taken 
into account for the photo-ionization rates but was found to 
be insignificant for the present case. The main effect of the 
line blanketing on the statistical equilibrium of metals that 
do not contribute significantly to the blanketing is generally 
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through the reduction of the available photoionizing UV pho- 
tons, but with such a large ionization energy, oxygen is insensi- 
tive to changes in the UV solar flux (see, e.g. Allende Prieto 
et al. 2003). Flux profiles were constructed from the angle- 
dependent intensity profiles. The corresponding LTE-profiles 
for the selected snapshots were also produced with MULTI3D 
using the same model atom but with all collisional probabilities 
multiplied by very large factors to ensure full thermalization of 
all rates. The non-LTE abundance corrections were estimated 
by computing both the non-LTE and LTE cases for three abun- 
dances differing by 0.2 dex and interpolating the non-LTE line 
strengths reproducing the LTE equivalent widths from the pro- 
file fitting. Both 3D snapshots gave the same non-LTE correc- 
tion to within 0.01 dex. Similarly robust results were obtained 
in the case of Li (Asplund et al. 2003a). Hence, while the line 
strengths vary somewhat from snapshot to snapshot following 
differences in the granulation pattern, the non-LTE abundance 
corrections stay essentially the same. 

The O model atom consists of 22 O I levels and the ground 
state of O II. The term diagram of the adopted model atom is 
presented in Fig. |2l In total 43 bound-bound and 22 bound- 
free radiative transitions were included together with elec- 
tron collisional excitation and ionization. Charge transfer reac- 
tions from the ground states have been taken into account; 
large charge transfer cross-sections for excited states as for Li I 
(Barklem et al. 2003) are not expected for O. We note that re- 
moving all charge transfer reactions only implies < 0.005 dex 
changes to the computed non-LTE abundance corrections for 
the Sun, reflecting the very minor influence of all ionization 
rates. Inelastic collisions with H for excitation and ionization 
were not accounted for, as the available laboratory measure- 
ments and detailed calculations strongly suggests that the clas- 
sical Drawin (1968) formula over-estimates the collisional rates 
by at least three orders of magnitude in the cases of Li and Na 
(Fleck et al. 1991; Belyaev et al. 1999; Barklem et al. 2003; see 
also discussion in Kiselman 2001). Unfortunately the specific 
case of O has not yet been investigated but there are no com- 
pelUng reasons on atomic physics grounds to expect a quali- 
tatively different behaviour for O compared with Li and Na in 
this respect. Furthermore, center-to-limb variation observations 
of the O I triplet suggests that the Drawin recipe over-estimates 
the collisional efficiency (Allende Prieto et al., in preparation). 
The adopted model atom is similar to the atom employed in 
the study of Kiselman (1993). We have verified that even an 
equivalent two-level atom yields very similar non-LTE abun- 
dance corrections as the full O atom. With the expected main 
non-LTE effect being photon losses in the line itself for the IR 
triplet our adopted atom is extensive enough to produce reliable 
results (Kiselman 1993). The 3D non-LTE abundance correc- 
tions for the O I lines employed in the present study are listed 
in Tabled 

Due to their high excitation potential (Xcxc > 9 eV), the O I 
lines are formed in deep atmospheric layers with a significant 
temperature contrast between the warm upflowing gas (gran- 
ules) and the cool downflows (intergranular lanes), as exem- 
plified in Fig. |3l It is clear that the main effect of the non- 
LTE line formation is a general strengthening of the Ol IR 
triplet over the whole granulation pattern. In this sense, the 



Table 2. The computed non-LTE abundance corrections for the 
three different model atmospheres. The departures from LTE 
have been estimated using MULTI3D (Botnen 1997 Asplund et 
al. 2003a) for the 3D case and with MULTI (Carlsson 1986) for 
the ID model atmospheres, in all cases using the same 23-level 
model atom. 
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-0.24 


-0.27 
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777.53 


5 go 
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9.15 


+0.00 


-0.20 


-0.24 


-0.20 


844.67 
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9.52 


+0.01 


-0.20 


-0.25 


-0.21 


926.60 


5p_ 


5po 


10.74 


+0.82 


-0.08 


-0.11 


-0.08 



non-LTE line formation is largely a ID problem: quite simi- 
lar results are obtained for O l if the non-LTE calculations are 
restricted to 1.5D (i.e. ignoring all non- vertical radiative trans- 
fer) rather than 3D. The depression of the line source function 
thus operates qualitatively similarly in both up- and downflows, 
although the details differ somewhat. As seen in Fig. |3l the 
largest 3D non-LTE effects are in fact seen in the intergranular 
lanes, which however, also house the smallest non-LTE effects 
in regions immediately adjacent to edge-brightened granules. 
Fig. |4] shows the line strength of the Ol 777.41 nm line as a 
function of the local continuum intensity. Again, the non-LTE 
and LTE behaviour are qualitatively similar for the triplet with 
the main difference being a shift in overall equivalent widths 
for the same abundance. The lack of a clear difference between 
the two cases in this respect is in sharp contrast to the case of 
Li (Kiselman 1997: Uitenbroek 1998; Asplund et al. 2003a), 
where the shapes of the LTE and non-LTE loci differ markedly. 
It will therefore be difficult to confirm the 3D non-LTE line 
formation modelling through observations of spatially resolved 
O I lines, as is possible for elements like Li and Fe (Kiselman 
& Asplund 2001). As expected, the [O l] lines are formed per- 
fectly in LTE, as seen in Fig.|4] 

The ID non-LTE line formation with the MARCS and the 
Holweger-Miiller model atmospheres has been studied with the 
code MULTI (Carlsson 1986) using the same model atom as in 
the 3D case. The structures of the MULTI3D and MULTI codes 
are very similar with the exceptions of the geometry and the use 
of long characteristics for the solution of the radiative transfer 
with MULTI. In particular, the subroutines for the background 
opacities and equation-of-state calculations are identical in the 
two cases. The non-LTE abundance corrections for the two ID 
model atmospheres are presented in Table|2] It is worth noting 
the similarity between the non-LTE corrections in ID and 3D in 
the case of O I lines in the Sun, in particular for the theoretical 
MARCS model. Although detailed results require 3D calcula- 
tions, non-LTE computations in ID model atmospheres are of 
great value to explore whether a given spectral line is prone 
to departures from LTE, at least at solar metalhcities. At low 
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Fig. 3. The granulation pattern in one of the solar snapshots 
seen in disk-center continuum intensity at 777 nm {upper left 
panel) and equivalent width of the O i 777.41 nm line in LTE 
{upper right panel) and non-LTE {lower left panel); the equiv- 
alent width images have the same relative intensity scale to 
emphasize the overall difference in line strengths. Also shown 
is the ratio of the non-LTE and LTE equivalent widths {lower 
right panel). 

metallicities, the temperature structure is sufficiently different 
between ID and 3D model atmospheres that one can expect 
significantly larger non-LTE effects in 3D for many elements 
(Asplund et al. 1999a, 2003a). 

3. The solar photospheric O abundance 

3. 1. Forbidden [Oi] lines 

The ability of the 3D hydrodynamical solar model atmosphere 
to very accurately predict detailed spectral line shapes and 
asymmetries without invoking additional ad-hoc broadening 
like micro- and macroturbulence represents a major advan- 
tage in abundance analyses, as obvious from the analysis of 
the forbidden [Ol] 630.0 nm line. It has long been suspected 
that the line is indeed blended by a weak Ni I line (Lambert 
1968). Due to the until recently very uncertain transition prob- 
ability of the Ni line, it has been difficult to ascertain its over- 
all contribution to the 630.0 nm feature and hence the solar O 
abundance. Lambert (1978) estimated the Ni contribution to 
be 0.02 pm=0. 2 mA, while the use of the theoretical 5f/-value 
(log 5/ = -1.74) of Kurucz (1993) for the Ni line would im- 
ply a very small left-over contribution for the [O l] line using 
the solar Ni abundance derived from other isolated lines. The 
3D model atmosphere can be used to estimate the necessary 
Ni line contribution in order to reproduce the shape of the ob- 
served feature. Without taking into account the Ni line, the pre- 
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Fig. 4. Predicted intensity (/i = 1.0) equivalent widths of the 
[Ol] 630.03 nm {left panel) and Ol 777.41 nm {right panel) 
lines across the granulation pattern in LTE {upper panel) and 
non-LTE {middle panel). Both cases are computed for log eo = 
8.90 for the two lines. The mean intensity equivalent widths 
are denoted by horizontal lines. Also shown are the ratios of 
the non-LTE and LTE equivalent widths {lower panel), with 
the horizontal lines representing the average ratio for the two 
lines. As expected, the [Ol] line is in perfect LTE while the 
O I line shows significant departures from LTE. Note that for 
clarity only a selection of the vertical atmospheric columns are 
shown. 



dieted [O i] fine has completely the wrong central wavelength 
and a line shape different from observed (see Fig. 1 in AUende 
Prieto et al. 2001), a discrepancy which is far worse than seen 
in similar 3D calculations for any other unblended line. 

From a x^-analysis of the flux profile of the [O i] 630.0 nm 
line, Allende Prieto et al. (2001) concluded that the solar O 
abundance is log eo = 8.69 ± 0.05 based on the same 3D 
model atmosphere as employed here. Since then, Johansson et 
al. (2003) have measured a new 5/- value for the blending Ni i 
line and details of its isotopic spfitting. Johansson et al. have 
modelled the shape of the Ni i line at the lab with two major 
components, due to the isotopes ^°Ni and ''^Ni, which con- 
tribute 94.3 % of the terrestrial abundance (Rosman & Taylor 
1998). Although the exact 51 /-value is inconsequential for our 
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purposes (only the product eNi • gf enters the analysis^), the 
additional line broadening from isotopic splitting changes the 
line profile of the Ni line somewhat and therefore also slightly 
the x^-analysis. The center of gravity of the line is also mildly 
changed from the value used by AUende Prieto et al. (630.0339 
nm; Litzen et al. 1993) to 630.0341 nm. Taking this into ac- 
count, however, only results in a 0.005 dex higher abundance 
than without the isotopic splitting. We therefore still arrive at 
logeo = 8.69 ± 0.05 (Table |3j- The same abundance is ob- 
tained when using disk-center intensity profiles instead of flux 
profiles. 

Adopting the theoretical line strength of only the [Ol] 
630.0 nm line as judged by the 3D analysis (0.36 pm=3. 6 mA 
for disk-center intensity and 0.43 pm for flux) for the ID cal- 
culations using the Holweger-Miiller and the MARCS model 
atmospheres, we find logeo — 8.77 ± 0.05 and logeo = 
8.73 ± 0.05, respectively. It should be noted, however, that 
it is not possible to accurately evaluate the Ni line contribu- 
tion to the feature based only on a ID analysis (Lambert 1968; 
Reetz 1998) due to the absence of intrinsic line asymmetries 
and the poor overall agreement between predicted and observed 
line profiles with ID hydrostatic models. For a long time, it 
was therefore simply assumed that the Ni line was unimpor- 
tant, leading to a much larger derived solar O abundance (e.g. 
Lambert 1978; Sauval et al. 1984): logeo - 8.90 for the 
Holweger-Miiller model using the new g/-value of the [O l] 
630.0 nm line adopted here. This clearly illustrates the impor- 
tance of having a realistic model atmosphere which adequately 
describes the intrinsic line shapes and asymmetries due to the 
convective motion and temperature inhomogeneities. 

As mentioned above, the analysis of the [Ol] 636.3 nm 
line is hampered by blending CN lines. In addition, the line 
is located in the midst of a Ca I auto-ionization line. The main 
blending CN line is the Q2 (25.5) 10-5 line with a center wave- 
length of 636.3776 nm, which is very close to the wavelength 
of the [O I] line: 636.3792 nm. In addition, the Pi (53.5) 8-3 CN 
line at 636.3846 nm is blending in the red wing. The contribu- 
tions of these CN lines have been determined from spectrum 
synthesis using the parameters for the two bands determined 
from other CN lines as well as for theoretical values. Together 
the two CN lines are estimated to have an equivalent width 
of 0.050 ± 0.005pm (flux), where the uncertainty reflects the 
different possible choices of transition probabilities of the CN 
lines. This leaves 0.14pm (1.4 mA) to be attributed to the [O l] 
line. In terms of abundance this corresponds to log eo = 8.67 
when analysed with the 3D hydrodynamical solar model at- 
mosphere (Table|3}. For the Holweger-Miiller and the MARCS 
model atmospheres, we find log eo = 8.75 and log eo — 8.71, 
respectively. 

The [Ol] 557.73 nm line is not considered here since it is 
badly blended by a C2 doublet whose exact contribution is very 
difficult to ascertain (Lambert 1978). 



^ In principle, one should re-determine the solar Ni abundance us- 
ing other Ni lines with our new 3D model atmosphere in order to eval- 
uate whether the new g/-value for the Ni I 630.0 nm line is indeed 
consistent with the solar Ni abundance. We postpone such a study to 
a later time as it is not crucial for our conclusions here. 



3.2. Permitted Oi lines 

Of the 20 odd observable lines attributed to O I in the wave- 
length range 500 — 1400 nm, unfortunately only a few are 
suitable for accurate abundance analysis besides the two for- 
bidden lines discussed above. The Ol 615.68 nm line is def- 
initely blended as are five lines at wavelengths greater than 
900nm (926.27, 974.15, 1130.24, 1316.39 and 1316.49nm). 
The heavy blending preventing their use has been verified by 
detailed 3D spectrum synthesis, as the shapes of clean lines 
are very well predicted by the new generation of 3D hydrody- 
namical model atmosphere (e.g. Asplund et al. 2000b; Asplund 
2000; Allende Prieto et al. 2002a). The two very faint (Wx < 

0. 1pm) lines at 645.44 and 700.19nm are not considered here 
because of their weakness and their uncertain theoretical gf- 
values. The 844.63 nm line is affected by uncertainties in the 
continuum placement and blending while the 976.07 nm lacks 
appropriate atomic data. This leaves only the familiar O I IR 
triplet at 777 nm (777.19, 777.41, 777.53 nm) and the even 
more high-excitation 615.81 nm line. In addition, we include 
the 844.67 nm line but with half weight compared with the 
other lines as it is partly blended and located in a region with 
a less well-determined continuum. The latter line also shows 
similarly large non-LTE abundance corrections as the IR triplet, 

1. e. ~ —0.2 dex. Finally, we include the 926.60 nm line, also 
with half weight as the line is blended in the red wing by a tel- 
luric H2O line. As a consequence, it is not possible to use the 
flux atlas to derive the abundance. While the Brault & Neckel 
(1987) disk-center intensity solar atlas^ still shows a prominent 
telluric feature, the Liege Jungfraujoch (Delbouille et al. 1973) 
disk-center atlas^ has a diminished problem with telluric ab- 
sorption in this wavelength region. For this line only we there- 
fore fit the disk center intensity profile instead of the flux pro- 
file for this particular line. Note that the line consists in fact 
of three components at 926.5827, 926.5927 and 926.6007 nm 
with log (7/ = —0.72, +0.13 and +0.71, respectively, which 
have been taken into account for the theoretical calculations. 

The derived O l-based abundances have been obtained from 
the temporally and spatially averaged 3D LTE flux profiles and 
taking into account the departures from LTE calculated in Sect. 
12.51 We note that almost identical abundances are derived when 
relying on disk-center intensity profiles instead. The 3D LTE 
profiles have been computed for 100 snapshots corresponding 
to 50min solar time for in total 17 different rays. Through a 
comparison with the observed profiles, the individual abun- 
dances best reproducing the line shapes and strengths for the 
different Ol lines were estimated. To these 3D LTE abun- 
dances, the 3D non-LTE abundance corrections presented in 
TableEl are added. The final 3D non-LTE abundances are found 
in Table |3l Due to the in general poorer agreement between 
predicted and observed profiles with ID model atmospheres, 
the same profile fitting technique is not as easily employed in 
the ID cases. Instead, the theoretical equivalent widths in 3D 
computed with the thus obtained abundances have been used 
as "observed" equivalent widths to be reproduced with the two 



^ ftp.hs.uni-hamburg.de/pub/outgoing/FTS-atlas 
' http;//bass2000.obspm.fr/solar_spect.php 
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Table 3. The derived solar oxygen abundance as indicated by 
forbidden [O l] and permitted O I lines. The results for the O I 
lines include the non-LTE abundance corrections presented in 
Table |2] which have been computed specifically for the three 
different model atmospheres 
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The abundances derived using the 3D model atmosphere have 
been obtained from profile fitting of the observed lines. 
^ The predicted flux equivalent widths (except for 926.60 nm which 
is for disk center intensity) with the 3D model atmosphere using 
the best fit abundances shown in the fourth column. 
The derived abundances with the Holweger-Miiller (1974) and the 
MARCS (Asplund et al. 1997) ID hydrostatic model atmospheres 
in order to reproduce the equivalent widths presented in the fifth 
column. Note that for the [O I] 630.03 nm line use of ID models 
prevents accurate estimation of the significant Ni I line blending 
(Allende Prieto et al. 2001). Without the 3D results at hand, a 
significantly larger O abundance of log eo — 8.9 would therefore 
be obtained. 

Half weight compared with the other lines is assigned to the 
844.67 nm line due to problems with blending and continuum 
placement. 

" Half weight compared with the other lines is assigned to the 
926.60 nm line due to problems with blending and continuum 
placement. Note that the line in fact consists of three components, 
which has been taken into account in the spectrum synthesis. In 
order to minimize the effect of the telluric H2O blend, the disk- 
center intensity profile has been fitted, and hence the tabulated 
equivalent width is for the intensity rather than for flux profile as 
for all other lines here. 



different ID model atmospheres. This ensures that the same 
line strengths are obtained in both ID and 3D and isolates the 
impact of the new generation of 3D solar model atmosphere 
compared with existing ID models. As in 3D, the ID non-LTE 
abundance corrections individually computed for the two ID 
model atmospheres have been added to the results presented in 
TableEl 

It is noteworthy that the 3D non-LTE abundances in Table|3l 
show very gratifying agreement in spite of the very significant 
differences in non-LTE abundance coiTections between the dif- 
ferent lines. In particular, the 615.81 nm and 926.60 nm lines 
show a small departure from LTE (— O.OSdex and — O.OSdex, 
respectively, in terms of abundance) while the other four O I 
lines have substantial non-LTE effects (—0.20... — 0.27 dex), 
yet the final non-LTE abundances agree to within 0.06 dex. 
While the non-LTE effect on the 615.81 nm line is small, the 
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Fig. 5. Upper panel: The temporally and spatially averaged 
LTE flux profile from the 3D model atmosphere (+) together 
with the observed profile (solid line) for the 1 777.53 nm line. 
The theoretical profile has been computed with log eo = 8.86 
and convolved with a Gaussian coiTesponding to the instrumen- 
tal resolution. Clearly the 3D LTE profile lacks effects of de- 
partures from LTE in the line formation. Middle panel: The 3D 
non-LTE (solid line) and the 3D LTE (dashed line) line pro- 
files computed for one snapshot with MULTI3D (Botnen 1997 
Asplund et al. 2003a). The 3D non-LTE profile has been com- 
puted with log eo = 8.70 while for the LTE case an abundance 
of logeo = 8.90 was adopted. In terms of equivalent widths 
the two profiles have the same line strengths. Lower panel: 
The temporally and spatially averaged 3D profile taking into 
account departures from LTE (+) shown together with the ob- 
served profile (solid line). The 3D non-LTE profile has been 
obtained by multiplying the 3D LTE profile in the upper panel 
with the ratio of non-LTE and LTE profiles presented in the 
middle panel. While not a perfect substitute for the temporally 
averaged 3D non-LTE profile it nevertheless shows a very en- 
couraging agreement with observations. The observed line pro- 
files have been corrected for the solar gravitational redshift of 
633ms~^. 

3D LTE effect (3D LTE - ID LTE) is significantly larger than 
for the other O I lines in view of its very high excitation poten- 
tial (Xcxc = 10.74 eV) and large atmospheric formation depth. 
In sharp contrast, with the Holweger-Miiller model atmosphere 
the 615.81 nm line implies a significantly larger abundance 
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than the 777 nm triplet by ^ 0.15 dex. With the MARCS model 
atmosphere the agreement is better but still significantly poorer 
than in 3D. It should be noted however that the here performed 
ID analyses benefit greatly from the 3D analysis. The equiva- 
lent widths estimated from the 3D line profile fitting are much 
more accurate than those directly measured from the observed 
profiles using for example Gaussian or Voigt functions as com- 
monly done due to the asymmetries of the lines and the pres- 
ence of blends. We interpret the excellent agreement in 3D as 
yet additional evidence for a very high degree of realism in the 
3D model atmosphere and 3D line formation modelling. 

Not only are the derived abundances very consistent in 3D, 
the predicted profiles also closely resemble the observed pro- 
files, as evident from Fig. |5] While the 3D LTE profiles are 
broader and have a too pronounced blue asymmetry compared 
with what the observations indicate, including the 3D non-LTE 
effects results in almost perfect matches with the observed pro- 
files for the O I IR triplet. It should be noted that due to com- 
putational constraints the 3D non-LTE calculations can only be 
performed for a small number of snapshots, preventing proper 
temporally averaged 3D non-LTE profiles as can be straight- 
forwardly achieved in 3D LTE. Instead the effects of depar- 
tures from LTE on the 3D profiles have been included by mul- 
tiplying the temporally averaged 3D LTE profiles with the ra- 
tios of the 3D non-LTE and 3D LTE profiles computed with 
MULTI3d for the two selected snapshots. This procedure is 
justified since the differences in 3D non-LTE abundance cor- 
rections between different snapshots are expected to be small 
(Asplund et al. 2003a), as also confirmed by specific calcula- 
tions for O. Thereby the final 3D non-LTE profiles incorpo- 
rates the effects of the Doppler shifts introduced by convec- 
tive motions and solar oscillations, which in ID is assumed 
to be replaced with a depth-independent and isotropic micro- 
and macroturbulence. The agreement between predicted and 
observed profiles is in general significantly poorer in ID even 
when including an individually optimized macroturbulence pa- 
rameter. The excellent agreement between observed and pre- 
dicted profiles for the 3D non-LTE calculations is a very strong 
support for our modelling. 

Further support for the 3D model atmosphere and the 3D 
non-LTE line formation modelling is obtained from a compar- 
ison of the center-to-limb variation of the O I 777 nm triplet. 
As discussed in detail in Kiselman & Nordlund (1995), the ob- 
served dependence of line strength on disk location /i is much 
more shallow than predicted when assuming LTE regardless 
of whether a 3D or ID model atmosphere is employed. In 
contrast, the non-LTE case reproduces the observed behaviour 
very nicely, at least when not including the very large cross- 
sections for inelastic collisions with H implied by the classi- 
cal Drawin (1968) formula (see discussion in Kiselman 2001). 
Not surprisingly given the similarity in the employed 3D solar 
model atmospheres and the robustness of the non-LTE effects 
as discussed in Sect. 12.51 we confirm the findings of Kiselman 
& Nordlund (1995). In this respect, our much more extended 
model atom compared with their two-level approach is of little 
consequence, as the dominant non-LTE effect is a depression 
of the line source function due to escape of 777 nm line pho- 
tons (Kiselman 1993). As shown in Fig. |6l the 3D non-LTE 
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Fig. 6. The observed center-to-limb variations of the equivalent 
widths of the Ol 777 nm triplet are shown as bands with es- 
timated internal errors. The solid lines denote the average of 
two snapshots of the theoretical 3D non-LTE results for the 
three lines, while the dashed lines represent the corresponding 
3D LTE case. The here shown 3D LTE curves are very simi- 
lar to the average of all 100 snapshots computed with the same 
3D LTE line formation code used for the calculations of the 
[O l] and OH lines. As two snapshots are insufficient to yield 
an accurate abundance estimate, the 3D non-LTE results have 
been interpolated to logeo ~ 8.7 to overlap with the obser- 
vational estimates to high-light the excellent agreement in the 
center-to-limb variation in 3D non-LTE. In all cases, the LTE 
curves correspond to a 0.2 dex higher abundance than the non- 
LTE cases. Clearly, LTE fails to reproduce the observational 
evidence, while the non-LTE calculations are in good agree- 
ment. 

results for the center-to-limb variations of the triplet agree very 
well with the observational evidence, in particular for ^ > 0.5. 
We note that an even more impressive agreement for fi < 0.5 
would have been achieved had we relied on the only, to our 
knowledge at least, previously published center-to-limb vari- 
ation observations of the triplet by Miiller et al. (1968) and 
Altrock (1968). Our new solar observations are in excellent ac- 
cordance with those of Miiller et al. (1968) and Altrock (1968), 
except for the data points at /i = 0.17, the only disk position 
with ^ < 0.5 we acquired spectra for. 

3.3. OH vibration-rotation lines 

The availability of high-resolution solar atlases in the IR 
with exquisite signal-to-noise ratio such as the here employed 
ATMOS atlas enables the use of the OH vibration-rotation 
around 3 /im and the OH pure rotation Unes beyond 9 /im as 
oxygen abundance indicators (Grevesse et al. 1984; Sauval et 
al. 1984). The high-quality observed spectrum combined with 
the relatively small amount of problematic blending in the rele- 
vant wavelength regions ensures that highly precise equivalent 
widths (initially in milli-Kayser = milli-cm^^ due to the na- 
ture of the solar atlas and subsequently converted to pm) of the 
OH lines can be measured. For the analysis of the OH X^Tl 
vibration-rotation lines we utilise 70 fines in the P^-branch of 



M. Asplund et al.: Solar line formation: IV. The photospheric O abundance 



11 



the (1,0) and (2,1) bands with TV" < 20; the Q and R Hnes fall 
in wavelength regions completely dominated by CO2 and H2O. 
The (1,0) and (2,1) lines are all weak and hence insensitive to 
the non-thermal broadening (convective motions in 3D and mi- 
croturbulence in ID). Due to the large number of lines, ourOH- 
based oxygen abundances are determined using the measured 
equivalent widths with the assumption of LTE in the line for- 
mation. For vibration-rotation lines like these, LTE is likely an 
excellent approximation (Hinkle & Lambert 1975). The issue 
of non-equilibrium molecule formation for OH has garnered 
relatively little attention in the literature but the available evi- 
dence suggests that it is not a serious worry for the OH lines in 
the Sun at least (Sanchez Almeida et al. 2001; Asensio Ramos 
et al. 2003). 

The resulting O abundances with the 3D model atmosphere 
show encouragingly small scatter and without any trends with 
neither wavelength, excitation potential nor line strength, as 
clear from Fig. Considering the high degree of temperature 
sensitivity of OH molecule formation and the minute fraction 
of oxygen tied up in OH, this excellent agreement strongly 
suggests that the temperature structure in the 3D model at- 
mosphere is very close to reality in the relevant atmospheric 
depth interval (typical mean optical depth of line formation 
lo 

§''"500 ~ —1.4 in the case of the Holweger-Miiller model. 
The concept of typical line formation depth is of course less 
well-defined in the case of 3D models). Giving equal weights 
to all 70 hnes, the average of the individual abundances yields 
logeo =8.61 ±0.03. 

A similar analysis based on the ID Holweger-Miiller model 
atmosphere gives a significantly larger abundance: log eo — 
8.87 ± 0.03. While the scatter is only marginally larger than 
in 3D, there are disturbing trends in oxygen abundances as 
function of wavelength and excitation potential (Fig.|8}; these 
trends are clearly inter-related as the low excitation lines have 
shorter wavelengths. It should also be noted that there is a 
significant trend with line strengths when using equivalent 
widths measured in wavenumber units rather than wavelengths. 
The relatively strong trend with excitation potential is particu- 
larly revealing, implying that the temperature structure of the 
Holweger-Miiller model is inappropriate. In fact, the same con- 
clusion led Grevesse & Sauval (1998) to experiment with a 
modified temperature structure in order to remove the exist- 
ing trends but thereby also affecting the mean OH-based oxy- 
gen abundance. In this respect, the ID MARCS model atmo- 
sphere is performing better. The resulting trends are barely sig- 
nificant (Fig.|9j and the scatter is very similar to the 3D case. 
The MARCS-based abundance is log eo = 8.74 ± 0.03. 

3.4. OH pure rotation lines 

The pure rotational lines of OH X^H were first identified in 
the solar atlas by Goldman et al. (1981, w = 0), Blatherwick 
et al. (1982, w = 1 and w = 2) and Sauval et al. (1984, 
V = 3). Goldman et al. (1983) and Sauval et al. (1984) have 
advocated for their use as a prime abundance indicator of the 
solar oxygen abundance given the large number of clean lines 
beyond 9/Ltm. Our analysis utilises 127 of the weakest OH pure 
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Fig. 7. The derived solar oxygen abundance (filled circles) from 
OH vibration-rotation lines using the 3D hydrodynamical time- 
dependent simulation of the solar atmosphere (Asplund et al. 
2000b) as a function of wavelength, lower level excitation po- 
tential and line strength (in pm). The solid hnes denote least- 
square-fits giving equal weights to all lines. 

rotation lines belonging to the v = 0, v = 1, v — 2, and 
V = 3 states with N" = 23 — 43, whose equivalent widths 
have been measured from the ATMOS solar IR atlas. Care has 
been exercised to only consider OH lines minimally affected 
by telluric absorption. While the equivalent widths and line 
depths for the pure rotation lines are in general smaller than for 
the vibration-rotation lines, the pure rotation lines are partly 
saturated since they are formed in significantly higher atmo- 
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Fig. 8. Same as Fig.Qbut using the Holweger-Miiller (1974) 
semi-empirical model atmosphere and a microturbulence of 
Cturb = 1.0 km s^^ (filled circles with solid lines denoting 
least-square-fits). The open circles and dashed lines correspond 
to a larger microturbulence (^turb = 1.5kms~^) 



Fig. 9. Same as Fig. 0but using a theoretical MARCS model 
atmosphere (Asplund et al. 1997) and a microturbulence of 
Ctuib ~ 1.0 km s^^ (filled circles with solid lines denoting 
least-square-fits). The open circles and dashed lines correspond 
to a larger microturbulence (^tuib = 1.5kms~^) 



spheric layers (typical mean optical depth of line formation 
—2.0 ^ logTsoo ^ —1.6 in the case of the Holweger-Miiller 
model). Hence the pure rotation lines are sensitive to the con- 
vective Doppler shifts (supposedly represented by microturbu- 
lence in ID analyses) in addition to the temperature sensitivity 
normal for all molecules. 

The equivalent width-based O abundances derived with the 
3D model atmosphere show a small scatter, in particular for 
the hnes with W\ = 1.5 — 10 pm. Due to observational prob- 



lems in accurately measuring the line strengths of the weak- 
est lines, the scatter increases significantly for W\ < 1.5 pm. 
However, there are distinct trends with excitation potential and 
line strength for lines with Xcxc < 2eV and W\ > 10 pm, 
as seen in Fig. ^| These two trends are directly related since 
the lowest excitation lines are also the strongest. It is expected 
that the realism of the 3D model atmosphere starts to deteri- 
orate towards higher atmospheric layers, in particular due to 
the employed radiative transfer treatment (total radiative heat- 



M. Asplund et al.: Solar line formation: IV. The photospheric O abundance 



13 



ing/cooling determined by relatively few lines, lack of non- 
LTE effects, no consideration of Doppler shifts for the energy 
balance etc). These lines are also sensitive to the atmospheric 
velocity field. Less weight should therefore be placed on the 
stronger rotational lines. In the following, the analysis is re- 
stricted to the pure rotation lines with Wa — 1-5 — 10pm. 
Giving equal weight to these 69 lines results in log eo — 
8.65 ±0.02. 

We have also performed corresponding calculations with 
the MARCS and Holweger-Miiller ID model atmospheres (Figs. 
[T2I and II 1> . The ID analysis using the MARCS model at- 
mosphere produces even more spectacular trends with exci- 
tation potential and line strengths with almost 0.2 dex higher 
abundance for the strongest lines compared with the weakest 
lines (Fig. I12> . Adopting instead a high microturbulence of 
^turb = 1.5kms^^ instead of the standard ^turb = l.Okms"^ 
only goes about half-way in resolving the problems. Since 
such a high microturbulence is inconsistent with similar analy- 
ses of other species, we conclude that the temperature struc- 
ture in these high atmospheric layers in the MARCS model 
is clearly inappropriate for the analysis of OH pure rotation 
lines. The average abundance for the same 69 lines as for 
the 3D analysis is log eo — 8.83 ± 0.02. In this respect, the 
Holweger-Miiller (1974) model is performing remarkably well 
with essentially no trends with the standard microturbulence 
of i^turb — 1.0 km s^^ and a very small scatter (Fig. II It . 
This is somewhat surprising given the problems encountered 
for the OH vibration-rotation lines with the same model atmo- 
sphere (Sect. lXSt . However, the pure rotation lines are formed 
in higher layers than the vibration-rotation lines (typically 
—2.0 ^ logTsoo ^ —1.6 and logrsoo — —1.4, respectively, 
in the case of the Holweger-Miiller model). The average abun- 
dance with the Holweger-Miiller model is log eo = 8.82±0.01. 

As mentioned above, it is slightly surprising that the 
Holweger-Miiller model performs better than the 3D model for 
the pure rotation lines while the opposite is true for the vibra- 
tion lines. The main effects must reside in the differences in 
temperature structure and/or velocity field (convection in 3D 
and microturbulence in ID). The mean temperature stratifica- 
tions differ in the relevant atmospheric layers, as illustrated in 
Fig-E It should be noted that the continuum optical depth scale 
at 3 /im (typical wavelengths for vibration lines) is very similar 
to that of 500 nm, while the continuum optical depth unity at 
10 fim (typical for rotation lines) occurs at T500 ~ —1.3. This 
also causes a reversal of the granulation pattern when viewed 
at these long IR wavelengths: the bright (warm) material is lo- 
cated above the intergranular lanes and the dark (cool) mate- 
rial lays above the normal granules (see Fig. 1 in Kiselman & 
Nordlund 1995). Since the typical line formation depths dif- 
fer for the two types of OH lines, it is possible to remove the 
abundance trends present in Fig.fTOjwithout tampering with the 
excellent consistency achieved in Fig. by slightly decreas- 
ing the temperatures outside T500 ~ —1.6. According to test 
calculations, the necessary temperature modifications are rel- 
atively small and fully within the uncertainties of the 3D tem- 
perature structure given the inherent approximations for the en- 
ergy balance in these high atmospheric layers. The main prob- 
lem appears to be the relative fine-tuning of the temperature 



Table 4. The derived solar oxygen abundance as indicated by 
forbidden [O l], permitted O I lines, OH vibration-rotation line 
and OH pure rotation lines. The results for the O I lines include 
the non-LTE abundance coiTections presented in Table |2l The 
quoted uncertainties only reflect the line-to-line scatter for the 
different types of O diagnostics while no doubt the total error 
is dominated by systematic errors. 



lines 




logeo 




3D 


HM 


MARCS 


[Ol] 


8.68 ± 0.01 


8.76 ±0.02 


8.72 ±0.01 


Oi 


8.64 ± 0.02 


8.64 ±0.08 


8.72 ± 0.03 


OH vib-rot 


8.61 ± 0.03 


8.87 ±0.03 


8.74 ± 0.03 


OH rot 


8.65 ± 0.02 


8.82 ±0.01 


8.83 ± 0.03 



alterations necessary in order not to perturb the vibration lines. 
It is also possible that the resolution of the discrepancies can 
be found in underestimated temperature inhomogeneities in the 
OH line forming region with the present 3D model. A detailed 
comparison of observed and predicted absolute disk-center in- 
tensities is expected to be able to shed further light on this issue. 
While the strongest rotation lines are sensitive to the velocity 
field, the required changes make this a rather implausible solu- 
tion according to various test calculations. 

3.5. Summary 

As seen above and summarised in Table |3 the different oxy- 
gen diagnostics used for the present study yield very consis- 
tent results when analysed with the 3D hydrodynamical solar 
model atmosphere. Clearly the situation is much less satisfac- 
tory when relying on classical ID model atmospheres. We see 
no compelling reason to prefer or exclude any of the individ- 
ual abundance indicators as the results are in all cases quite 
convincing. The quoted errors in Table |3 only reflect the line- 
to-line scatter (standard deviation rather than standard devia- 
tion of the mean as occassionally used for solar analyses) for 
the different types of lines. Using these standard deviations, the 
weighted mean becomes logeo ~ 8.66 ± 0.02. However, the 
total error budget is most certainly dominated by systematic 
errors rather than statistical. Possible such sources include the 
atmosphere model (in particular the temperature structure), the 
continuous opacity (due perhaps to non-LTE effects on the H~ 
molecule) and non-LTE effects on the line and molecule for- 
mation. To evaluate systematic eiTors is notoriously difficult. 
An attempt of doing so can be achieved by a comparison of 
the abundances implied by the different indicators, which in the 
3D case suggests that the remaining systematic eiTors are under 
reasonable control given the very different types of line forma- 
tion processes involved and their widely different temperature 
sensitivities. A reasonable estimate would be to assign an un- 
certainty of ±0.05 dex on the derived solar oxygen abundance. 
We therefore arrive at our best estimate of the solar oxygen 
abundance as: 

logeo = 8.66 ±0.05. 
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4. The solar photospheric Ne and Ar 
abundances 

The solar abundances of neon and argon are also directly af- 
fected by the new solar abundance of oxygen. Since spectral 
lines of Ne and Ar do not show up in the solar photospheric 
spectrum, their solar abundances have to be deduced from coro- 
nal matter. Oxygen, neon and argon particles are present in 
the solar wind as well as in solar energetic particles (SEP). As 
they are high FIP (First lonisation Potential) elements they are 
not subject to fractionation affecting the low FIP elements in 
coronal matter From the observed accurate SEP ratios, Ne/O = 
0.152 ± 0.004 and Ar/O = (3.3 ± 0.2) • lO^^ (Reames 1999; 
these values agree with the somewhat less accurate solar wind 
data as well as with local galactic ratios derived from hot stars 
and nebulae) we deduce revised values of the solar abundances 
ofNeandAr: logeNc = 7.84 ±0.06 and log eai- = 6.18±0.08 
respectively. These values are again much smaller than the pre- 
viously recommended values of 8.08 and 6.40, respectively, in 
Grevesse & Sauval (1998). The quoted errors account both for 
the measurement uncertainty in the element ratios in SEP and 
the estimated systematic error in the solar O abundance. 

5. Comparison with previous studies 

The derived solar oxygen abundances from the different indi- 
cators presented in Sect. |3] are particularly noteworthy for two 
reasons: the significantly lower O abundance by as much as 
0.25 dex compared with previously commonly adopted values 
(Anders & Grevesse 1989; Grevesse & Sauval 1998), and the 
excellent agreement for all different features, in stark contrast 
to what is achieved with ID model atmospheres. In particular 
the latter fact is a very strong argument in favour of the new 
generation of 3D hydrodynamical model atmosphere and their 
application to accurate abundance analyses. 

In spite of their weakness, the forbidden [Ol] 630.0 and 
636.3 nm lines have often been advocated to be the superior 
O abundance indicators, in particular the former line (Lambert 
1978). As both lines originate from the ground state, their 
line formations are guaranteed to be exceedingly well ap- 
proximated by LTE as verified by our non-LTE calculations. 
Until recently, it was believed that the contribution from the 
blending Ni I 630.0 nm line was negligible, leading to an oxy- 
gen abundance of logeo — 8.9 with the Holweger-Miiller 
semi-empirical model atmosphere (Lambert 1978; Sauval et al. 
1984; Holweger 2001). Strong support for this high abundance 
has come from analyses of the OH vibration-rotation (Grevesse 
et al. 1984) and pure rotation lines (Sauval et al. 1984), which 
show an impressively small scatter between different OH lines 
using the Holweger-Miiller model. Together this led to the con- 
clusion that the solar O abundance is logeo = 8.93 ± 0.03 
(Anders & Grevesse 1989). Since then the OH-based abun- 
dance has gradually come down by 0. 1 dex in view of new 
preliminary analyses using a Holweger-Miiller model atmo- 
sphere with a slightly modified temperature stratification to re- 
move existing trends with excitation potential and line strength 
(Grevesse & Sauval 1998, 1999) and improved (^/-values and 
measured equivalent widths. In the meantime, it has become 



clear that the Ni I blend for the [O l] 630.0 nm line may not be 
insignificant after all. Reetz (1998) made a very comprehensive 
investigation using both flux and intensity solar atlases at dif- 
ferent disk locations in an attempt to constrain the Ni contribu- 
tion. While no unanimous conclusion could be reached, Reetz 
concluded that the Ni line could affect the results so much that 
the solar O abundances may be as low as log eq = 8.75 when 
analysed with the Holweger-Miiller model. A more precise es- 
timate of the Ni line contribution was obtained by Allende 
Prieto et al. (2001) by employing the same 3D hydrodynam- 
ical solar model atmosphere as used here. In terms of abun- 
dance it was found that properly taking the Ni line into ac- 
count will lower the O abundance by about 0.1 3 dex. In com- 
parison, the impact of moving from the classical Holweger- 
Miiller ID model atmosphere to a 3D model was less impor- 
tant in this case, 0.08 dex, but the two effects compound each 
other by going in the same direction. Together with a new gf- 
value (same as used here), Allende Prieto et al. (2001) obtained 
log eo = 8.69 ± 0.05. The work presented here builds on the 
success of this study but extends it greatly by also perform- 
ing similar calculations for the other O diagnostics. It is al- 
ways worrisome when a result rests entirely on only one spec- 
tral line, in particular when the implication in terms of derived 
abundance is as profound as in this case and the line itself is 
significantly blended. 

While both the [O l] and OH lines for long pointed towards 
a high solar O abundance, the permitted high-excitation Ol 
lines were causing severe problems, in particular the IR triplet 
at 777 nm. It is well-known that the triplet is susceptible to de- 
partures from LTE, which increases the line strength and thus 
decreases the derived abundance. The exact magnitude of this 
non-LTE effect depends on the recipe, if any, for excitation by 
inelastic collisions with H. The exact choice of model atmo- 
sphere and other details of the model atom play a smaller role in 
this respect. The essentials of the line formation process of the 
triplet lines are in fact very well described even with a two-level 
approach in as far as the main non-LTE effect is photon-losses 
in the line itself, causing the line source function to deviate 
from the Planck value (Kiselman 1993). A smaller contribu- 
tion from a slight change in the line opacity also modulates the 
result however. Without accounting for H-coUisions, the non- 
LTE abundance corrections for these lines are about —0.2 dex 
for the Sun, leading to a solar O abundance of log eq — 8.7 
(Kiselman 1993). Due to the large discrepancy in the non- 
LTE case with the [Ol] and OH results, many authors have 
argued that the non-LTE calculations are faulty. Indeed the 
LTE-abundance is close to those provided by the other diag- 
nostics, providing in some people's eyes a mean to calibrate 
the poorly known H-collisions (Tomkin et al. 1992; Takeda 
1994). By necessity, such a procedure requires the H-collisions 
to be very large, often even larger than the classical Drawin 
(1968) recipe for such processes, thus enforcing efficient ther- 
malization of the line and a result close to the LTE-expectation. 
It is important to realise, however, that the available labora- 
tory and theoretical calculations unambigously suggest that the 
Drawin recipe over-estimates the H-collisions by about three 
orders of magnitude in the sofar only studied cases of Li and 
Na (Fleck etal. 1991; Belyaevetal. 1999; Barklemetal. 2003). 
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If the same conclusions hold for O then the H-collisions are 
completely negligible in the solar case. Furthermore, center-to- 
limb variations of the predicted LTE -profiles (or those obtained 
from non-LTE calculations with Drawin-like H-collisions) dif- 
fer markedly from the observational evidence, as seen in Fig.|6] 
(see also Kiselman & Nordlund 1995). 

The root of the problem with ID analyses has therefore 
been the poor agreement between the O abundances estimated 
from the [O l] and OH lines on the one hand and the O I lines on 
the other; with the Ni-induced downward revision of the [O l]- 
based abundance this statement must be modulated however. 
In an important paper, Kiselman & Nordlund (1995) attempted 
to resolve this long-standing discrepancy through a pioneering 
3D analysis of the O I triplet, the [O l] 630.0 nm line and a few 
OH vibration-rotation and pure rotation lines. They used a sim- 
ilar but less sophisticated 3D model atmosphere compared to 
the one employed here and restricted the calculations to only 
two different snapshots. A noteworthy finding of their study 
was that it may be possible to reconcile the different abun- 
dance indicators in a proper 3D analysis since the [Ol] and 
OH lines became stronger than in the Holweger-Miiller model 
atmosphere while the opposite was true for the triplet as in- 
dicated by an equivalent two-level non-LTE calculation. Due 
to the restricted numerical resolution and temporal coverage, 
they were unable to perform an accurate abundance analyses 
for the lines, but they concluded that "the solar oxygen abun- 
dance could very well be less than 8.80". Our 3D model atmo- 
sphere covering some 100 snapshots over a 50min solar-time 
sequence with an improved temperature structure and exten- 
sion to greater heights, allows us to derive significantly more 
accurate abundances and finally settle the issue. 

Recently, Holweger (2001) has independently re-visited the 
solar O abundance using the [Ol] 630. Onm and Ol lines; he 
did not study molecular lines. Holweger used the old 5 /-value 
for the forbidden line and neglected the Ni-blend and conse- 
quently found a high oxygen abundance: logeo — 8.92 us- 
ing the Holweger-Miiller model. His seven permitted lines sug- 
gested a significantly lower value, which, together with non- 
LTE abundance corrections, yielded logeo = 8.73 ± 0.04. 
The difference relative to our value with the Holweger-Miiller 
model atmosphere can be traced to differences in selected 
lines, adopted equivalent widths and much smaller non-LTE 
corrections (for completeness it should also be noted that 
Holweger uses intensity rather than flux profiles). In particular, 
Holweger adopts large H collision cross-sections, leading to 
small non-LTE abundance corrections. As described above, the 
evidence is stacked against such efficient H collisions. In addi- 
tion, Holweger added "granulation corrections" amounting to 
—0.02 dex on average to the derived abundance, finally arriving 
at logeo = 8.74±0.08 when including the [Ol] 630.0nmUne. 
It should be noted, however, that these granulation corrections 
represent the difference in derived abundances between full 2D 
spectrum synthesis calculations and corresponding ID com- 
putations using the ID average of the 2D model atmosphere 
(Steffen & Holweger 2002). Thus, the abundance corrections 
are designed to only reflect the effect of temperature inhomo- 
geneities but not the difference in mean temperature structure 
with for example the Holweger-Miiller model atmosphere. As 



a consequence, these granulation corrections can differ signifi- 
cantly both in magnitude and sign from those estimated by our 
method, which also includes the effects of different mean tem- 
perature structures between the ID and 3D models. As our 3D 
solar model atmosphere has successfully passed a range of ob- 
servational tests (e.g. Stein & Nordlund 1998: Rosenthal et al. 
1999; Asplund et al. 2000b; Asensio Ramos et al. 2003), we 
are confident that our approach is fully justified. 



6. Conclusions 

We have presented a determination of the solar oxygen abun- 
dance from an analysis of [O l], O I, OH vibration-rotation and 
OH pure rotation lines by means of a realistic time-dependent 
3D, hydrodynamical model of the solar atmosphere. All oxy- 
gen indicators yield very consistent results, in marked contrast 
with studies employing classical ID model atmospheres. The 
here derived solar oxygen abundance of log eo = 8.66 ± 0.05 
(Table ^ is significantly smaller than previously thought: the 
revision amounts to almost a factor of two (0.27 dex) compared 
with the often quoted value of log eo ~ 8.93 from Anders & 
Grevesse (1989), and still much smaller than the value of 8.83 
recommended more recently by Grevesse & Sauval (1998). 
This low abundance is strongly supported by the excellent 
agreement between abundance indicators of very different tem- 
perature sensitivities and line formation depths, and between 
the observed and predicted line shapes and center-to-limb vari- 
ations. It should be noted that not the whole difference with pre- 
vious results is attributed to the employment of a 3D model at- 
mosphere over classical ID models since the adoption of more 
recent ^/-values, more realistic non-LTE procedures, better ob- 
servations and a proper accounting of blends also play an im- 
portant role in this respect. 

The new low solar oxygen abundance derived with the 
new generation of 3D hydrodynamical model atmospheres (e.g. 
Stein & Nordlund 1998; Asplund et al. 1999, 2000b; Asplund 
& Garcia Perez 2001) resolves the problematic high metal- 
licity of the Sun compared with the solar neighborhood sug- 
gested by previous studies. The oxygen abundance presented 
here (O/H— 460 ± 50 • 10"^) is in very good agreement 
with that measured in the local interstellar medium for real- 
istic gas-to-dust ratios (Meyer et al. 1998; Andre et al. 2003): 
OfH— 410 ± 60 • 10^^. In addition, the solar oxygen abun- 
dance is now similar to those obtained from studies of nearby 
B stars (Cunha & Lambert 1994; Kihan et al. 1994). The 
most recent collation suggests that the mean B star value is 
O/H = 350 ± 130 • 10"*^ (Sofia & Meyer 2001) but there is a 
significant scatter between the existing B star analyses. Given 
the uncertainties we do not consider the remaining differences 
with our new solar oxygen abundance as significant. Finally, 
our value is close to those of nearby young. Galactic disk F and 
G dwarfs: 0/H= 450 ± 160 • lO"*^ (Sofia & Meyer 2001). 
Interestingly, while the previous discrepancy most often has 
been blamed on flaws in the analyses of hot stars and nebu- 
lae, the most serious shortcoming has apparently rather been on 
the solar side. We believe that 3D line formation calculations 
of late-type stars such as those presented herein finally address 



16 



M. Asplund et al.: Solar line formation: IV. The photospheric O abundance 



this important issue, promising much more accurate abundance 
analyses for the future. 

There is, however, at least one area where the new solar 
O abundance together with the corresponding downward re- 
visions of the photospheric C, N and Ne abundances (Asplund 
2003, 2004; Asplund et al. 2004b, c) does not resolve a problem 
but rather severely aggrevates it. The accurate measurement of 
the solar oscillations, helioseismology, allows a unique glimpse 
of the structure of the solar interior (Christensen-Dalsgaard 
2002). The analysis can be posed as an inversion problem, en- 
abling the extraction of the variation of the sound speed with 
depth in the Sun. The predictions from standard models of 
the solar structure and evolution show a remarkable agreement 
with the observed sound speed from helioseismology to within 
Acs/cs = 4 • 10^3 (e.g. Christensen-Dalsgaard 2002). The 
existing inversions are normally based on the opacities com- 
puted with the old Anders & Grevesse (1989) solar chemical 
composition, or at least compositions with much larger values 
than our new results. With the new significant down-ward revi- 
sions of the solar C, N, O and Ne abundances, the overall solar 
metallicity in the photosphere decreases from Zqm — 0.0194 
(Anders & Grevesse 1989) to Z^cw = 0.0126 (mass fraction of 
elements heavier than He). This diminished metallicity prop- 
agates to a significant modification of the opacity because C, 
N, O and Ne are main contributors to the opacity in different 
layers below the convection zone, and the sound speed is there- 
fore directly affected. As a result, the excellent agreement be- 
tween the predicted and observed sound speed variation with 
depth deteriorates markedly (Pijpers 2002, private communi- 
cation; see also Boothroyd & Sackmann 2003). We note that 
the new solar chemical composition should also affect the de- 
termination of the solar He abundance in the outer convective 
zone from helioseismology inversion due to modifications of 
the standard solar model. Likewise, the initial He abundance 
in the protosolar cloud will likely be slightly smaller than the 
value Y — 0.27 (mass fraction of He) normally derived with 
existing solar evolution models based on the old metal abun- 
dances. 

At this stage it is unclear whether a resolution to this prob- 
lem can be found and if so where the solution may lay. In the 
meantime, we offer the tentative suggestion that the problem 
may be traced to the fact that the here obtained solar C, N and 
O abundances only refer to the photospheric values rather than 
the abundances found in the solar interior If the diffusion of 
metals, as incorporated into the stellar evolution models used 
for the solar inversion, is more efficient than currently thought, 
it is conceivable that the good correspondence in terms of the 
sound speed variation could be recovered. The metallicity in 
the interior would then be significantly higher than that of the 
convection zone and the photosphere. The excellent concor- 
dance between all the different C, N and O diagnostics achieved 
with the new generation of 3D hydrodynamical solar model at- 
mospheres in contrast with ID analyses certainly strongly ar- 
gues for that the metallicity of the solar photosphere is indeed 
significantly lower than previously thought. 
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Fig. 10. The derived solar oxygen abundance from OH pure 
rotational lines using a 3D hydrodynamical solar model at- 
mosphere as a function of wavelength, lower level excitation 
potential and line strength (in pm). The solid lines denote the 
least-square-fits. The rotational lines are more sensitive to the 
adopted microturbulence than the vibration-rotation lines since 
the stronger lines are partly saturated. The increased scatter for 
the weakest lines (which correspond to the highest excitation 
lines) is most likely due to increased observational difficulties 
in measuring the equivalent widths accurately 



Fig. 11. Same as Fig.^Jbut using the Holweger-Miiller model 
atmosphere. The filled circles correspond to a microturbulence 
of 'Cturb ~ l.Okms^^ while the open circles are for ^turb = 
l.Skms^^ with the solid and dashed lines, respectively, denot- 
ing the least-square-fits 
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Fig. 12. Same as Fig.llOlbut using a MARCS model atmosphere. 
The filled cuxles correspond to a microturbulence of ^turb = 
1.0 km s^^ while the open circles are for i^turb — 1.5 km s^^ 
with the solid and dashed lines, respectively, denoting the least- 
square-fits 



